Microstructure Effects for Casimir Forces in Chiral Metamaterials 
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We examine a recent prediction for the chirality-dependence of the Casimir force in chiral meta- 
materials by numerical computation of the forces between the exact microstructures, rather than 
homogeneous approximations. We compute the exact force for a chiral bent-cross pattern, as well as 
forces for an idealized "omega" -particle medium in the dilute approximation and identify the effects 
of structural inhomogeneity (i.e. proximity forces and anisotropy). We find that these microstructure 
effects dominate the force for separations where chirality was predicted to have a strong influence. 
To get observations of chirality free from microstructure effects, one must go to large separations 
where the effect of chirality is at most ^ 10~^ of the total force. 



It has been proposed that dielectric metamaterials 
might exhibit repulsive Casimir forces in vacuum where 
planar structures have only attraction [IHl]- However, 
these predictions used effective-medium approximations 
(EMAs) for the metamaterials, often treating the EMA 
terms as free parameters. While certain parameters 
give repulsion, a recent theorem [5 implies that one 
cannot have repulsion between vacuum-separated di- 
electric/metallic metamaterials in the effective-medium 
(large separation) regime (although it is possible outside 
of this regime [6]). However, it is possible that signif- 
icant reduction or modulation of the Casimir force can 
occur as a result of metamaterial effects. A recent EMA 
analysis of chiral metamaterials [7 predicted such a force 
reduction, and even repulsion, due to chirality. Unfortu- 
nately, chirality effects vanish at large separations where 
the EMA should be valid, and are strongest at small 
separations where the EMA is questionable. Therefore, 
the question of whether chiral metamaterial effects can 
theoretically modulate the Casimir force remains open. 
Answering it requires accurate calculations using the ex- 
act microstructure of the metamaterials, rather than a 
homogeneous EMA. Recent advances [8UT4] have made 
the treatment of such complex structures possible. In 
this paper, we apply these methods to rigorously test the 
EMA predictions for chiral metamaterials against numer- 
ical calculations incorporating the microstructures. We 
are able to distinguish chiral metamaterial effects from 
other "non-ideal" effects such as pairwise surface-surface 
attractions, and find that the former are overwhelmed 
except at large separations, where they are only 10~^ or 
less of the total force. 

Dielectric and metallic metamaterials are defined by an 
inhomogeneous permittivity £:(x, cj), but at sufficiently 
long wavelengths (large separations in the Casimir con- 
text) they can be accurately modeled in the EMA by 
homogeneous constitutive parameters [e.g, £{uj)^ 



that can be very different from the constituent materials. 
The Casimir force between two bodies, however, is nat- 
urally expressed as an integral over imaginary frequency 
uj = [15 . The force in the EMA will then depend on 
the effective e{i(,), /^(^O^ ^^c. over a range of imaginary 
frequencies [4 . Early works [1-3 predicted repulsive ef- 
fects by a putative €{i^) < /i(^0 body. However, 
no such repulsion was found for "magnetic" metamate- 
rials based on actual structures, in which e{i(,) > /i(iO 
for all ^ [4^ 16^. Basically, while these metamaterials can 
have almost any e and /i at a given real uj via resonances, 
on the imaginary-cj axis the important features come 
from the behavior at low ^ (long wavelengths). In this 
limit, /i(zO - m(0) ^ -^'^ gimilT], where /i(0) < ^(O). 
Therefore, there is no repulsion or force reduction in the 
EMA regime for this model. More recently, Ref. 7 stud- 
ied chiral metamaterials, which in addition to e and /i are 
characterized in the EMA by a chirality coupling D to 
H and B to E. ^ — by a spatial inversion (x — x) 
of the microstructure. For any the EMA predicts that 
the force between two media of the same chirality (SC) 
should be lower than for media of the opposite chirality 
(OC), with the size of this difference increasing with k,. 
Because in general ^ ^ ^Tj ilTj , this effect is largest at 
short separations and goes to zero at large separations. 
Therefore, to get a useful prediction (e.g., force reduction 
due to chirality) , we are forced to consider the predictions 
of the EMA at intermediate separations. However, not 
only does the EMA fail at sufficiently small separations, 
but the force in that limit is eventually described by the 
proximity force approximation (PFA) [18 in which there 
are only pairwise attractions. One should therefore be 
cautious of any EMA prediction that differs qualitatively 
from PFA in this limit. An analysis based on pairwise 
attractions would find little effect due to chirality. In- 
stead, the behavior would be dominated by small-scale 
features such as the shortest distance between two com- 
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FIG. 1: (color online) Interactions between two chiral meta- 
materials [19]. Shaded box indicates a unit cell; the structures 
are periodic in the transverse direction. Nearest neighbor bars 
on opposing structures are colored yellow. In the EMA, the 
forces should obey Fsc < Foe- By contrast, in the pairwise- 
force approximation the force is determined by the number of 
overlapping nearest-neighbor bars, as indicated: Foe > Fsc 
when the centers are aligned (left column), but Foe < Fsc 
when they are displaced by L/2 (right column). 

ponents of the microstructures. Thus far, no attempt has 
been made to determine which behavior dominates, and 
we do this as follows: if the materials behave as homoge- 
neous, chiral media, the relative chirality should be the 
only source of a force difference between the SC and OC 
cases, where Fqc > Fsd independent of the transverse 
displacement x. However, if the force is governed by pair- 
wise attraction, it should exhibit a strong x-dependence. 
We can therefore directly test the validity of the EMA 
by comparing the force for different values of x. 

We first examine a realistic structure proposed 
in Ref. and shown in Fig. [T] A unit cell of each 
"medium" consists of a single bi-layer of two bent-arm 
crosses, one spiraling clockwise and the other counter- 
clockwise. Their ordering in the 2:-direction determines 
the chirality of the medium. We omit the dielectric poly- 
imide in which the metal was embedded in Ref. 19 , as this 
eliminates a chirality-independent attraction between the 
layers. The exact Casimir force between the periodic 
structures of Fig. [l] is computed using a finite-difference 
time-domain (FDTD) method ^Q^. Ten bi-layers in the 
2;-direction are included on each side; adding more layers 
does not change the results. We compute the force for 
two different material types: perfect electric conductors 
(PEC) and dispersive gold. For gold we take a plasma 
model with ujp = 1.37 x lO^^rad/sec. The latter model 
requires a definite value for the length scale a, which we 



FIG. 2: (color online) Relative force difference {Foe — 
Fse) / Foe for the left (red) and right (blue) columns of Fig. IT] 
for PEC (solid) and dispersive gold for a = 1 /im (dashedjT^ 
From Fig. [l] the sign difference is due to proximity effects. 
The inset shows a zoom of the large- z regime, along with the 
EMA prediction (black). 

take to be a = 1/im. The results shown in Fig. |2] (each 
force difference is normalized by Foc{x)) are similar for 
both PEC and dispersive gold, and are not consistent 
with the EMA: even the sign of Foc{x) — Fsci^) can be 
changed as a function of x for z/a < 0.75. Larger z/a 
still exhibit a strong x-dependence in the force. 

These results for low z/a can be qualitatively described 
by pairwise nearest-neighbor attractions. Consider, for 
instance, the two cases diagrammed in Fig. [T] First, 
when the centers of the unit cells are aligned, there is 
approximately a 4-bar overlap in the SC case and an 8- 
bar overlap in the OC case. If the force is proportional 
to the number of overlapping nearest-neighbor bars, we 
expect Foe > Fsc- When the centers are displaced by 
X = L/2, SC has a 4-bar overlap and OC has a 3-bar 
overlap, so we expect Fqc < Fsc^ and the relative force 
difference should be reduced by approximately 1/4. This 
prediction, based purely on pairwise attraction, captures 
the behavior for z/a < 0.75. For larger z/a the sign of 
Fqc — Fsc is the same in both cases, and the magnitude 
of the force difference is comparable to the EMA predic- 
tion. However, in this limit chiral effects are very small 
(accounting for only 0.1% of the force). Furthermore, as 
the force difference is still highly x-dependent, and other 
transverse displacements may still switch the sign of the 
force difference. We therefore cannot determine if this is 
really an "ideal" chiral effect or not. 

As we have seen, when the chirality of the microstruc- 
ture is thought to have a large effect on the force, this 
effect can actually be attributed to pairwise forces rather 
than chiral metamaterial effects. However, the predic- 
tions of Ref. [19] were based on isotropic^ chiral media, 
and the bent-arm crosses of Fig. [l] are highly anisotropic. 
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FIG. 3: (color online) Top: Unit cell of isotropic chiral meta- 
material: PEC "omega" particles in vacuum. Shaded cube 
indicates relative positions in unit cell. Bottom: Effective- 
medium parameters, as a function of imaginary frequency 
cj = i^, deduced directly from imaginary-cj scattering data. 
/i(0) / 1 because the particles are PEC. 



It is possible that the effects of this anisotropy domi- 
nate those of chirahty. To obtain a chirahty with min- 
imal anisotropy, we now consider a more idealized sys- 
tem: the so-called "omega" particles [17, 19 , shown 
in Fig. [3| Each individual particle was predicted to 
have a strong chiral response, and when assembled into 
a period-a unit cell of high symmetry should form an 
isotropic chiral metamaterial. The unit cell configura- 
tion is shown in Fig. [s] (top). Fig. [s] (bottom) shows 
the obtained EM A parameters (discussed below), making 
this an isotropic chiral metamaterial. In fact, by picking 
a unit cell of lower symmetry (see below), we find that 
anisotropy quickly dominates chiral metamaterial effects. 
To make the medium as homogeneous as possible, the 
omega particles are relatively closely packed, and as dis- 
persion had little qualitative role in the previous results, 
we take the particles to be PEC for simplicity. 

Force computations for a structure as complex as Fig.jS] 
are very difficult for finite-difference methods, because of 
the disparity of size between the periodicity a and the 
wire diameter 0.016a. Instead, we employ a boundary- 
element formulation [10 that computes the scattering 
matrices of objects in imaginary uj using a nonuniform 
mesh of the surfaces, shown in Fig.[3j In addition, we use 
a dilute approximation (justified later) in which multiple 



scattering events within a given structure are neglected, 
so only the scattering matrices of individual particles are 
required. The force between two periodic structures is 
then computed from the scattering matrices (see Ref. JT] 
for a detailed derivation, and a partial review of pre- 
cursors p!2Hl4] ). The force computation is summarized 
as follows: first, the scattering matrix for each omega 
particle is numerically computed in a spherical multipole 
basis. This leads to a matrix of scattering amplitudes 
Fu ,m' ,P';i,m,p^ whcrc —I < TTi < I and —I' <m'<l' 
are the respective moments of the incident and scattered 
spherical waves, and P, P' G {M, E} their polarization. 
Second, the scattering matrix is converted to a planewave 
basis, where for each uj the planewave has transverse 
wavevector and polarization P = p. Third, to get 
the scattering from the unit cell of Fig.[3]we sum the scat- 
tering matrices from the twelve individual omega parti- 
cles, each rotated and displaced appropriately. The rota- 
tions are applied to the spherical multipole moments via 
a rotation of the spherical harmonics, and the displace- 
ments are applied in the planewave basis via planewave 
translation matrices. Fourth, the two-dimensional peri- 
odicity of the lattice is incorporated into the problem: 
for scattering from a unit cell in the planewave basis, one 
computes the scattering coefficients between planewaves 
of arbitrary and k'^, P and P^ When the scattered 
fields are summed over all unit cells of a periodic struc- 
ture (in two dimensions), the scattering matrix gains a 
factor Ec^^^Hkx - + G)/a^, where {G} are the 
reciprocal lattice vectors. Therefore, only scattering be- 
tween planewaves where k^ — k'^ = G are needed. Fi- 
nally, to simulate a semi-infinite omega medium in 2:, we 
apply the translation matrices to the scattering ampli- 
tude of an entire unit cell, displaced by integer multiples 
of a in With this method, we can quickly compute 
forces for many configurations, e.g., many x — z displace- 
ments. For the present computations, we find that / < 3 
and k^ within the first three Brillouin zones suffice to 
get the force (and the force difference) to high precision. 

To demonstrate that we have an isotropic, chiral 
medium in the EMA, we extract the effective-medium 
parameters s{i^)^ /^(^O^ ^(^0 [Fig. [s] (bottom)] from 
the scattering matrix at k^. Parameter retrieval [19] 
from reflection and transmission at normal incidence can- 
not be used, because for imaginary u the transmission de- 
creases exponentially with the thickness of the medium. 
Instead, we compute the specular reflection coefficients 
Rss{oJ^'k±) and Rsp{oJ^'k±). For each frequency uj = i^, 
these quantities (given in Ref. 7) can be expanded to 
quadratic order in <C ^ and tz. The coefficients of 
each term are determined by a fit (with |kx| < ^/lO). 

The results of the force computations are shown 
in Fig. [4] and are similar in form to Fig. [2] We exam- 
ine the relative force difference at each z/a — normalizing 
by the minimum force is convenient to obtain a positive 
difference that can be plotted on a log scale. The shaded 
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FIG. 4: (color online) Forces (in dilute approximation) be- 
tween chiral media from Fig.js] for the SC (blue) and OC (red) 
cases. The range of the force across all transverse (x) displace- 
ments is shaded, and normalized by the minimum of the force 
over all x. For z < 3.1a, chirality is not well-defined (curves 
overlap). Only when curves are distinct and x-independent 
can the systems be described as chiral metamaterials. In- 
set: Frequency integrand of the relative force difference at 
z — 3.6a. The smallness of the u — term is an indicator 
that this difference is due to chirality and not anisotropy (see 
text). 



areas of each color represent the range of values that the 
force assumes for all transverse {x) displacements. Their 
spread (the x-dependence) indicates breakdown of the 
EMA, and for z < 3.1a it is so severe that the red and 
blue force curves overlap. When these curves separate, 
but before they sharpen (the blue (SC) curve goes to 
zero), we have an intermediate regime where chiral ef- 
fects are competing with proximity effects. Only when 
the curve thickness is much less than the relative force 
difference (< 10% for z > 3.6a) is EMA accurate. As 
mentioned above, we designed the unit cell to be highly 
isotropic, for which the EMA parameters ji^n are truly 
scalars. For a less isotropic unit cell (e.g., only one omega 
particle per unit cell), the relative force difference can 
be made much larger, but the source of this difference 
is ambiguous, as in Fig. [2] For two purely chiral ma- 
terials in the EMA regime, the zero- frequency compo- 
nent of the force difference is exactly zero [7], so the 
magnitude of the zero-frequency component is an indica- 
tion of whether the force difference arises from chirality 
alone. In Fig. [Ij^Inset) we display the frequency inte- 
grand at z = 3.6a, indicating that the force difference 
in the isotropic case is indeed due to chirality. How- 
ever, we have found that the frequency integrand for less 
isotropic unit cells has a large zero- frequency component, 
indicating that anisotropy and proximity effects play a 
large role in that case. This could also explain the re- 
sults of Fig. |2] for larger z. In conclusion, to observe an 
unambiguous effect of chirality in isotropic, chiral meta- 



materials, one must measure the total force to 4 digits 
of accuracy. Given that a 1 fim for realistic systems, 
and that anisotropics in the structure can result in much 
stronger force variations, this effect is unlikely to be rel- 
evant for applications. This does not preclude the pos- 
sibility of other interesting Casimir effects in metamate- 
rials (especially if they rely on the u ^ response of 
the metamaterial), but it imposes serious limitations for 
finite-o; effects. 

Finally, we comment on the validity of the effective 
medium picture and the dilute approximation. One can 
compute the Casimir force and the force difference from 
the EMA parameters of Fig. [s] For the range of z/a 
in Fig. [4| this gives the correct force between the two 
media but overestimates the force difference by roughly 
30 percent. This is because the error terms in the EMA 
above are 0{kj_), the same order as the chirality contribu- 
tions ~ i?gp to the force [7 . The Clausius-Mossotti (C- 
M) equation allows us to check the validity of the dilute 
approximation. C-M relates polarizabilities of particles 
to the effective- medium parameters ^ . We com- 
pute the polarizabilities from C-M, working backwards 
from the effective-medium parameters of Fig. [s] (expand- 
ing C-M to first order in the polarizability) , and plug 
this into the full C-M to compute the non-dilute correc- 
tion. The relative force difference of Fig. [4] is changed 
by < 15%. C-M is obtained in the static limit, whereas 
at finite imaginary u the exponential field decay reduces 
the interactions, so this is an overestimate. 
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